Hot electrons in low-dimensional phonon systems 
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A simple bulk model of electron-phonon coupling in metals has been surprisingly successful in ex- 
plaining experiments on metal films that actually involve surface- or other low-dimensional phonons. 
However, by an exact application of this standard model to a semi-infinite substrate with a free sur- 
face, making use of the actual vibrational modes of the substrate, we show that such agreement 
is fortuitous, and that the model actually predicts a low-temperature crossover from the familiar 
T 5 temperature dependence to a stronger T 6 log T scaling. Comparison with existing experiments 
suggests a widespread breakdown of the standard model of electron-phonon thermalization in metals. 

PACS numbers: 63.22.+m, 85.85.+j 
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The coupling between electrons and phonons plays 
a crucial role in determining the thermal properties of 
nanostructures. The widely used "standard" model of 
low temperature electron-phonon thermal coupling and 
hot-electron effects in bulk metals assumes (i) a 
clean three-dimensional free-electron gas with a spheri- 
cal Fermi surface, rapidly equilibrated to a temperature 
T e i; (ii) a continuum description of the acoustic phonons, 
which have a temperature T p h; (iii) a negligible Kapitza- 
like thermal boundary resistance |3| between the metal 
and any surrounding dielectric, an assumption that is of- 
ten well justified experimentally; and (iv), a deformation- 
potential electron-phonon coupling, expected to be the 
dominant interaction at long- wavelengths. In a bulk 
metal, the net rate P of thermal energy transfer between 
the electron and phonon subsystems is 



p = Ev e in 



T, 



where V c \ is the volume of the metal, and 
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Here £ is the Riemann zeta function, cf is the Fermi en- 
ergy, N e \ is the electronic density of states (DOS) per unit 
volume, p is the mass density, v\ is the bulk longitudinal 
sound speed, and vf is the Fermi velocity. 

This model, which has no adjustable parameters, has 
successfully explained some experiments 0, 0] but 
others report a power-law temperature dependence with 
smaller exponents |6|,l3> indicating an enhanced electron- 
phonon coupling at low temperatures. However, the 
experiments typically involve heating measurements in 
thin metal films deposited on semiconducting or insulat- 
ing substrates, and the relevant phonons at low temper- 
ature are strongly modified by the exposed stress-free 
surface. An attempt to directly probe such phonon- 
dimensionality effects was carried out by DiTusa et at 
0, who intentionally suspended some of their samples, 
necessarily modifying the vibrational spectrum, although 



they found no significant difference from their supported 
films. We argue that the paradox reported in Ref. 0] is 
actually quite widespread, and all experiments known to 
us on supported films actually contradict the standard 
model when that model is modified to account for the 
actual vibrational modes present in a realistic supported- 
film geometry, illustrated in Fig. ^ Our results have im- 
portant implications for the thermal properties of meso- 
scopic and low-dimensional phonon systems and the use 
of such systems as nanoscale thermometers, bolometers, 
and calorimeters 0, . 

The Hamiltonian we consider (suppressing spin) is 

H = J2k £ k C k c k+Z)„ fr^n ala n +5H, where and c k are 
electron creation and annihilation operators, with k the 
momentum, and and a n are bosonic phonon creation 
and annihilation operators. The vibrational modes, la- 
beled by an index n, are eigenfunctions of the continuum 
elasticity equation v% V x V x u — vfV( V • u) = lu 2 u for 
linear isotropic media, along with accompanying bound- 
ary conditions. Vt and V\ are the bulk transverse and 
longitudinal sound velocities. 5H = |ep fydPriptyV ■ u 
is the deformation-potential electron-phonon interaction, 
with u(r) = J]„(2pwn)~^[fn(r)an + f*(r)a£] the quan- 
tized displacement field. The vibrational eigenfunctions 
f„(r) are defined to be solutions of the elasticity field 
equations, normalized over the phonon volume V^h ac- 
cording to f v d?r f* • f„/ = 8 nn t. It will be conve- 
nient to rewrite the electron-phonon interaction as SH = 
Ekqn [ffnq 4+qCk d„ + s£ q 4_ q c k 4] , with coupling con- 
stant g nci = le F (2puj n )-iV c { 1 J v J 3 r V • f„ e - iC i r . Note 
that we allow for different electron and phonon volumes. 

The quantity we calculate is the thermal energy per 
unit time transferred from the electrons to the phonons, 



P ee 2 ]T Hu n [ir (k - k - q) - r» b (k • 

kqn 
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where 



rT(k - k - q) = 2t: \g nci \ 2 [n B (w„) + 1] 
x n F (e k )[l - np(ek-q)] <S(ek- q - £k + w n ) 



(4) 
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is the golden-rule rate for an electron of momentum k to 
scatter to k — q while emitting a phonon n, and 

r* b (k -> k + q) = 2tt |g, iq | 2 n B K) 

x riF(e k ) [1 - n F (ek+q)] #(ek+q - £k - w n ) (5) 

is the corresponding phonon absorption rate, ne is the 
Bose distribution function with temperature T p h and np 
is the Fermi distribution with temperature T e \. The fac- 
tor of 2 in (|3J) accounts for spin degeneracy. It is possible 
to obtain an exact expression for P; the result (suppress- 
ing factors of % and fes) is 
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The logarithmic term in P can be shown to be negli- 
gible in the temperature regime of interest and will be 
dropped. Carrying out the k integration then leads to 



P = 



vf£V el 



duj F(uj) 



24C(5bo 

where F(w) = J2 n U n <5( w 
brational DOS, with 
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Here wd is the Debye frequency. U n can be interpreted 
as an energy associated with mass-densit y fl uctuations 
interacting via an inverse-square potential 11], cut off at 
distances of the order of the lattice constant a. We have 
reduced the calculation of P to the calculation of F(uj). 
Allen ^3 nas derived a related weighted-DOS formalism. 

We now calculate F(oj) and P for a metal film of thick- 
ness d attached to the free surface of an isotropic elastic 
continuum with L — > oo; see the inset to Fig. ^ F° r 
calculational simplicity, the film and substrate are as- 
sumed to have the same elastic parameters, characterized 
by a mass density p and bulk sound velocities Vt and v\. 
Where material parameters are necessary we shall assume 
a Cu film; however, the qualitative behavior we obtain is 
generic. The evaluation of F(u>) requires the vibrational 
eigenfunctions for a semi-infinite substrate with a free 
surface, which have been obtained in the classic paper 
by Ezawa [uj- The modes are labeled by a branch in- 
dex m, taking the five values SH, +, — , 0, and R, by 
a two-dimensional wavevector K in the plane defined by 
the surface, and by a parameter c with the dimensions 
of velocity that is continuous for all branches except the 
Rayleigh branch m = R. With the normalization conven- 
tion of Ref. 0| we have 



F(lu) = %k S(lo - c R K) 

dc UmKc S(U! - cK). 
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FIG. 1: (inset) Conducting film of thickness d attached to 
insulator. The top surface of the metal is stress-free, (main) 
Temperature dependence of the thermal power exponent x for 
a 10 nm (solid curve) and 100 nm (dashed curve) Cu film. 



The range of the parameter c depends on the branch m, 
and is summarized in Table [I] The frequency of mode 
mKc is cK. 



m 


range of c 


SH 


[v t , oo] 


± 


[v\, oo] 





[vt , V\] 


R 


cr (discrete) 



TABLE I: Values of the parameter c for the five branches of 
vibrational modes of a semi-infinite substrate. 

Turning to an evaluation of (JSJl, the SH branch is purely 
transverse, so J7sh = 0. The normalized eigenmodes for 
the ± branches are 



f ± = 



AttcA 1 



=F ol 2 (e 



C±e' 



x ( e - ll3Kz + C±e ll3Kz ) 



ex 



± q5 ( e ~* aKz + f ± e laKz ) 



+ i/ri(e-^*-C±e^*) 
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where a = y/ (c/v\) 2 — 1 and /3 = \J(c/v%) 2 — 1. Here 
[(/3 2 - 1) ± 2i^a~p] 2 



{(3 2 - l) 2 + 4a/3 



with |c±|=i. (10) 



Then 
V • f± = Ti 
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and U± = {(?K/av?V e i)I±(Kd,c), where 



and U R = {2-Kd^K/K.vfVex) I R (Kd), where 



r z , r z 

I±(Z,c) = RcJ dxdx' K ( y /(x-x') 2 + a 2 Z 2 /d 2 ) J R (Z) = J dxdx' K (y/(x-x') 2 + a 2 Z 2 /d 2 )e-^ x+x '\ 



x r e ia(x-x') _ A_gia(x+x')l 



Ko is a modified Bcssel function. To obtain J7± we use 
translational invariancc in the xy plane to write J7J as 



U m Kc = 77 / dzdz' [ d 2 R 
*4i Jo J A 



i? 2 + (z - z') 2 + a 2 ' 

where R = (x, y) is a two-dimensional coordinate vector. 
Then we scale out K, do the angular integration, and use 
the identity /" dRRl (R) [R 2 + s 2 ]^ 1 = Ko(M), where 
J is a Bcssel function of the first kind. 

Next we consider the m — branch, for which 



iC e-~< Kz + i(3 e- lf3Kz + i(3A e lf3Kl 



ex 



7 C e~^ Kz + ie~ tl3Kz ~ 1A e l(5Kz 



; }e* Kr , (13) 



where 7 = y/l— (c/v{) 2 , 
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Then 



V-f = 



e ' e 



and J7 = (\C\ 2 c 3 K/(3vfV cl )I (Kd,c), where 



(14) 



I (Z,c) = dx dx' K (^(x-x') 2 + a 2 Z 2 jd 2 ) e-^+x") . 
Jo 

Finally, for the Rayleigh branch, 
f R = 



K 
KA 



-ipKz 



Spy \ „-v Kz 



( 2yr; \ 
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where <p = a/1 - (c R /vi) 2 , 77 = ^/l - (c R /i> t ) 2 , and /C = 
(if — T))((p — 77 + 2tp7f)/2ipi] 2 . cr is the velocity of the 
Rayleigh surface waves, given by cr = £vt, where £ is the 
root between and 1 of £ 6 8£ 4 + 8 (3 - 2 ^ 2 ) £ 2 - 1 6 ( 1 - z; 2 ) , 
with v = Vt/v\. For Cu, v — 0.52 and £ = 0.93; hence 
cr = 2.4 x 10 5 cnis _1 . Using fT5|) . 



'£4 



(16) 



The final summations in (jSJ) are carried out with the 
aid of the identity lim^^oo 5{u> — cK) = u>A/2irc 2 
and elsewhere replacing K with lu/c. Then we obtain 

+ £°° <ic ii j+( ^'' ?)+ - f - ( ^' c) ]}- (i7) 

This expression, combined with ©, is our principal re- 
sult. Evaluation of (|17[) can be further simplified by the 
use of the powerful identities 



2? 



' df(Z.s) 
\ ds 



I ± (Z,c) = Re[(2Z-^)f(Z,a) + ^e 2 * aZ r(Z.a) 

(18) 

lf(Z,-i 7 ), (19) 
-f(Z-iip), (20) 



I (Z,c) = I/(Z,i 7 )-; 
7r(Z,c) = lf{Z,iy)- 



where /(Z, s) = $*dxK Q (y/x* + a 2 Z 2 /d 2 ) e lsx , thereby 
reducing the /,„ to a single one-dimensional integral /. 

The I m have distinct large- and small-Z character, 
crossing over near Z = 1 . Because of the integration over 
c in l|17fl . F and P accordingly exhibit a broad crossover 
behavior. However, once ojd < cr, all branches will have 
assumed their low-frequency forms. We define a crossover 
temperature 



T* = hc K /k B d 



(21) 



dividing regimes determined by the small and large 
tod/cYi behavior of F. In the large wc?/cr limit the 
m = ± modes in (|17J) can be shown to be dominant, and 
lim W( j^co f^dc ^/±(^,c) = wuid- Therefore, we obtain 
F(u) — > Fbuik(^) = ^ 3 /wf, independent of c?, leading to 
a high-temperature behavior P — > S V^i [$(cjD/Iei) T e 5 j — 
$(wD/7ph) , where S is the coefficient (J2J, and where 

$(y) ee [4! ((SjrVo' dx ^VO* _ X )- $ ( 10 ) is about °- 97 ' 
and $(y) rapidly approaches 1 beyond that. Thus, at 
temperatures above T* but sufficiently smaller than the 
Debye temperature, the $ factors are equal to unity, and 
we recover the bulk result QJ. 

The low temperature asymptotic analysis is somewhat 
complicated and will be presented elsewhere. Briefly, us- 
ing the small Z expansion 



f{Z,s) -> -Z\nZ+ (l + ln2 + V(l))Z- —Z 2 \nZ 
+ ^(i+ln2 + ^(l))Z 2 + 0(Z 3 lnZ), (22) 
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where ip is the Euler polygamma function, we find 
F{lo) -> F bulk (o;)x[-A(^)ln(^) + 0(^)] in the small 
Lod/cn limit. Here 



2ttc 2 /3 



dc 



c R [2-Re(C++^)] 
2nc 2 a 



is a constant determined by v\, Vt, and cr,. Each T 5 func- 
tion in Q therefore crosses over at low temperature as 
T 5 -A(^) In (^r) , with A = A tt 6 /189 C(5). For a Cu 
film, A « 0.815 and A sa 3.998. There are also mixed- 
temperature regimes possible, where only one of the two 
terms in has crossed over. 

The most striking consequence of the crossover is that 
the temperature exponent increases. In Fig. ^ wc fit P 
(with either T B \ or T p h zero) to a power-law T x with a 
temperature dependent exponent x, and plot the expo- 
nent for 10 nm (T* = 1.84 K) and 100 nm (T* = 184 mK) 
Cu films. x(T) is nonmonatonic, displaying a pronounced 
maximum near T*, and drifts upward as T — > 0. Such 
behavior has not (to the best of our knowledge) been ob- 
served, even though many experiments 0, [1 EL have 
achieved T <C T*. The physical origin of the crossover 
is that, at low temperature, the stress-free condition at 
the metal surface penetrates into the film, reducing the 
strain and hence electron-phonon coupling there. The 
characteristic distance over which the boundary condi- 
tion has an effect is of the order of a bulk wavelength. 
When T>T*, only a thin outer surface layer of the film 
has a significantly diminished strain, and bulk behavior 
is expected. However, when T <C T* the entire metal 
film experiences a reduced strain. 

The experiments of Refs. 4] and both using 
Cu films, observe an approximate T 5 dependence even 
well below T*. It is therefore interesting to compare 
the observed pref actors with the coefficient S, evalu- 
ated for Cu. Using a free-electron gas ap proximation 
[lij and measured elastic properties |l5l |. we obtain 
5.97 x 10 7 W m -3 K -5 , which is at least an order of mag- 
nitude smaller than observed, consistent with our asser- 
tion that there is some unidentified mechanism enhanc- 
ing the thermal coupling. Noble metals are far from 
free-electron systems because of their complex Fermi sur- 
faces. We attempt to address this shortcoming by re- 
garding the "Fermi surface" quantities N c i(e-p) and vf 
as independently adjustable parameters, to be obtained 
empirically from heat capacity and cyclotron resonance 
data. Carrying out this analysis, the details of which 
will be presented elsewhere, leads to the modified prefac- 
tor £ = 1.14xl0 8 Wm~ 3 K~ 5 , which is still considerably 
smaller than measured. 

Although not included in the model considered here, 
disorder in a bulk metal film is expected to produce a 
crossover from the T 5 dependence to a T 6 scaling when 
the phonon elastic mean free path £ becomes smaller than 
the thermal wavelength 0, 0], a behavior which has 



not been reported experimentally until very recently . 
Thus, the crossover predicted here should not be appre- 
ciable affected by disorder unless £ < d. Although thin 
films are known to scatter phonons strongly, measured 
values of I are still much larger than d in the tempera- 
ture regime of interest here [19j. 

In conclusion, we argue that a wide variety of exper- 
iments contradict the predictions of an essentially ex- 
act application of the standard model of electron-phonon 
thermal coupling in metals to a supported-film geometry, 
suggesting a widespread breakdown of that model. ANC 
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and CMS-040403. 



[5. 
[6 

[-; 

[8 
[9 

[10 

[11 

[12 
[13 
[14 



V. F. Gantmakher, Rep. Prog. Phys. 37, 317 (1974). 
F. C. Wellstood, C. Urbina, and J. Clarke, Phys. Rev. B 
49, 5942 (1994). 

W. A. Little, Can. J. Phys. 37, 334 (1959). 
M. L. Roukes, M. R. Freeman, R. S. Germain, R. C. 
Richardson, and M. B. Ketchen, Phys. Rev. Lett. 55, 
422 (1985). 

C. S. Yung, D. R. Schmidt, and A. N. Cleland, Appl. 
Phys. Lett. 81, 31 (2002). 

J. Liu and N. Giordano, Phys. Rev. B 43, 3928 (1991). 
J. F. DiTusa, K. Lin, M. Park, M. S. Isaacson, and J. M. 
Parpia, Phys. Rev. Lett. 68, 1156 (1992). 
M. L. Roukes, Physica B 263, 1 (1999). 
K. Schwab, E. A. Henriksen, J. M. Worlock, and M. L. 
Roukes, Nature (London) 404, 974 (2000). 

D. R. Schmidt, R. J. Schoelkopf, and A. N. Cleland, 
Phys. Rev. Lett. 93, 45901 (2004). 

Recall that, in elasticity theory, the mass-density fluctu- 
ation Sp is given by —pV ■ u. 
P. B. Allen, Phys. Rev. Lett. 59, 1460 (1987). 
H. Ezawa, Ann'. Phys. (N.Y.) 67, 438 (1971). 
In the free-electron gas approximation, the Fermi energy 
of Cu is 7.03 eV, the DOS at e F is 1.13 x 10 34 erg" 1 cm" 3 , 
and «f is 1.57 x 10 s cms -1 . 
[15] Cu has the low-temperature elastic constants cu = 

1.25 x 10 12 dyne cm" 2 , and 



1.76 x lO^dynecm- 3 , c 12 = 
C44 = 8.18 x lO^dynecm" 



|2C|. and its mass density 
is p = 8.94 gem -3 . To approximate a cubic material as 
an isotropic elastic continuum we define effective elastic 
constants Cu = cu — 2A and C44 = C44 + A, where 
A = (en — C12 — 2c44)/5, and we let vi = (Cu/ p)? and 
ft = {Cn/p)^. This leads to V\ = 4.97 x 10 s cms -1 . 
[16] B. Keck and A. Schmid, J. Low. Temp. Phys. 24, 611 
(1976). 

[17] J. Rammer and A. Schmid, Phys. Rev. B 34, 1352 (1986). 
[18] I. J. Maasilta, J. T. Karvonen, J. M. Kivioja, and L. J. 

Taskinen, e-print cond-mat/0311031. 
[19] T. Klitsner and R. O. Pohl, Phys. Rev. B 36, 6551 (1987). 
[20] W. C. Overton and J. Gaffney, Phys. Rev. 98, 969 (1955). 



